Generalized virial theorem in f(R) gravity 
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We generalize the virial theorem in f(R) modified gravity using the collisionless Boltzmann equa- 
tion. We find supplementary geometric terms in the modified Einstein equation providing an effective 
contribution to the gravitational energy. The total virial mass is proportional to the effective mass 
associated with the new geometrical term, which may account for the well-known virial theorem 
mass discrepancy in clusters of galaxies. The model predicts that the geometric mass and its effects 
extend beyond the virial radius of the clusters. We also consider the behavior of the galaxy cluster 
velocity dispersion in f(R) models. The metric inside the galactic cluster, as well as the Lagrangian 
of the modified gravity model, are obtained in terms of quantities directly related to the physical 
properties of the clusters, and which can be determined from astrophysical observations. Thus, the 
f(R) virial theorem can be an efficient tool in observationally testing the viability of this class of 
generalized gravity models. 

PACS numbers: 04.50.+h, 04.20.Jb, 04.20.Cv, 95.35.+d 



(N 
> 

O 



o 



X 



I. INTRODUCTION 

The issue of dark matter is a long outstanding prob- 
lem in modern astrophysics. Two observational aspects, 
namely, the behavior of the galactic rotation curves and 
the mass discrepancy in clusters of galaxies led to the 
necessity of considering the existence of dark matter at 
a galactic and extra-galactic scale. The rotation curves 
of spiral galaxies show that the rotational velocities in- 
crease from the center of the galaxy and then attain an 
approximately constant value, vt goo ~ 200 - 300km/s, 
within a distance r from the center of the galaxy [l[ . In 
these regions the mass increases linearly with the radius, 
even where very little luminous matter can be detected. 
Relatively to the mass discrepancy in clusters of galaxies, 
the total mass of a cluster can be estimated in two ways. 
First, by taking into account the motions of its member 
galaxies, the virial theorem provides an estimate, My. 
Second, the total baryonic mass M may be estimated by 
considering the total sum of each individual member's 
mass. The mass discrepancy arises as one generally veri- 
fies that My is considerably greater than M, with typical 
values of M v /M ~ 20 - 30 Q. 

This is usually explained by postulating the existence 
of a dark matter, assumed to be a cold pressure-less 
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medium distributed in a spherical halo around the galax- 
ies. A number of candidates for dark matter have been 
proposed in the literature, in particular, the most pop- 
ular ones being the weakly interacting massive particles 
(WIMP) 0. Their interaction cross section with nor- 
mal baryonic matter, while extremely small, is expected 
to be non-zero and therefore may be directly detectable. 
However, despite more than 20 years of intense experi- 
mental and observational research, the non- gravitational 
evidence for dark matter is still lacking. Moreover, ac- 
celerator and reactor experiments do not yet support the 
physics beyond the standard model, on which the dark 
matter hypothesis is based upon. It is interesting to note 
that dark matter consisting of WIMPs may exist in the 
form of an Einstein cluster [3j , or could possibly undergo 
a phase transition to form a Bose-Einstein condensate [3] • 
It has also been suggested that the dark matter in the 
Universe might be composed of super-heavy particles, 
with mass > 10 10 GeV. But observational results show 
that dark matter can be composed of super-heavy parti- 
cles only if these interact weakly with normal matter or 
if their mass is above 10 15 GeV Q. 

Therefore, it seems that the possibility that Einstein's 
(and Newton's) gravity breaks down at the scale of galax- 
ies cannot be excluded a priori, and several theoret- 
ical models, have been proposed in the literature [f3]. 
One possibility is to assume that at large scales gen- 
eral relativity breaks down, and a more general action 
describes the gravitational field. The theoretical mod- 
els in which the standard Einstein-Hilbert action can be 
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replaced by an arbitrary function of the Ricci scalar R 
have recently been extensively investigated [7[. In this 
context, cosmic acceleration can arise due to corrections 
to the usual gravitational action of general relativity in 
the form R + 1/R, as now the 1/R term dominates for 
low curvatures |8|. Therefore, modified gravity, in prin- 
ciple, eliminates the need for dark energy, and since a 
modification of the Einstein-Hilbert action of the form 
f(R) oc R + R 2 can lead to early-time inflation Q, it 
should also be possible to unify the early and the late 
time accelerating phases of the Universe. 

However, it was shown that in all f(R) theories that 
behave as a power of R, at large or small R, the scale fac- 
tor during the matter phase grows as t x l 2 instead of the 
standard law i 2 / 3 This behavior, which is also a gen- 
eral property of models with d 2 f/dR 2 < 0, is inconsistent 
with cosmological observations (e.g. WMAP), thereby 
ruling out these models, even if they pass the supernovae 
test and can escape the local gravity constraints. Nev- 
ertheless, ,f(R) modified theories of gravity can, in gen- 
eral, give rise to cosmological viable models compatible 
with a matter-dominated epoch evolving into a late ac- 
celerated phase [ll|, [ijj. Thus, the viability of the f{R) 
models proposed has been a fundamental issue exten- 
sively analyzed in the literature [H, 03, EI El, Ell • 
In this context, severe weak field constraints in the so- 
lar system range seem to rule out most of the models 



Il8l . [2fJ , although viable models do 



proposed so far 

exist [Silil 

Specific models of f(R) gravity that admit a modified 
Schwarzschild-de Sitter metric and describe both the Pio- 
neer anomaly and the flat rotation curves of spiral galax- 
ies were analyzed in [24L l25j . In particular, within the 
framework of f(R) gravity, a model has been proposed 
exhibiting an explicit coupling of an arbitrary function of 
R with the matter Lagrangian density [261 ] , which estab- 
lishes a connection between the problem of the rotation 
curve of galaxies, via a solution somewhat similar to the 
one put forward in the context of MOND, and the Pio- 
neer anomaly. The possibility that the galactic dynamics 
of massive test particles may be understood without the 
need for dark matter was also considered in the frame- 
work of f(R) gravity [13, [H, [29[. The results obtained in 
these papers seem to suggest that a strong modification 
of standard general relativity is required to explain the 
observed behavior of the galactic rotation curves. How- 
ever, it was shown that to explain the motion of test 
particles around galaxies, in the constant rotational ve- 
locity regions, only mild deviations from classical general 
relativity are required [301 ] . 

Due to its generality and wide range of applications, 
the virial theorem plays an important role in astro- 
physics. Assuming steady state, one of the important 
results which can be obtained with the use of the virial 
theorem is to deduce the mean density of astrophysical 
objects such as galaxies, clusters and super clusters, by 
observing the velocities of test particles rotating around 
them. Hence the virial theorem can be used to predict the 



total mass of the clusters of galaxies. The virial theorem 
is also a powerful tool for stability studies. In a general 
relativistic framework several versions of the virial the- 
orem have been proposed 1311 , including the effect of a 
cosmological constant |32l |33j , and the generalization to 
brane world models [34j |. 

It is the purpose of the present paper to consider the 
virial theorem in the framework of f(R) modified theories 
of gravity. Using the collisionless Boltzmann equation in 
the modified Einstein field equations we derive a gener- 
alized virial equality. It takes into account the presence 
of the supplementary geometric effects due to the modi- 
fication of the gravitational action. 

This supplementary geometric terms give an effective 
contribution to the gravitational energy, the total virial 
mass being proportional to the effective mass associated 
to the new geometrical terms. Therefore this term may 
account for the well-known virial theorem mass discrep- 
ancy in clusters of galaxies. By using the virial theo- 
rem together with the gravitational field equations of the 
f(R) model the metric inside the cluster of galaxies can 
be obtained in a simple form, in terms of physical param- 
eters that can be fully obtained from observations, like, 
for example, the temperature of the intra-cluster gas and 
the radius and central density of the cluster core. Once 
the metric is known, the Lagrangian of the modified f(R) 
gravity model can also be obtained in terms of observable 
physical quantities. To explain the observational data at 
the extra-galactic scale a logarithmic modification of the 
standard general relativistic action is needed. Therefore 
the generalized virial theorem in f(R) gravity can be an 
efficient tool in observationally testing the viability of 
this class of generalized gravity models. 

This paper is outlined in the following manner. In 
Section ITT] we present the gravitational field equations of 
the galactic clusters in the scalar-tensor version of f(R) 
modified theories of gravity. The relativistic Boltzmann 
equation is outlined and the respective generalized virial 
theorem in f(R) gravity is deduced in Section HTT1 Astro- 
physical applications are explored in Section [IV! in par- 
ticular, predictions of the geometric mass and geometric 
radius from galactic cluster observations are presented, 
and the behavior of the galaxy cluster velocity disper- 
sion in f(R) models is also considered. In Section |V] we 
derive the Lagrangian of the modified gravity model, and 
analyze some of its properties. Finally, in Section PVTl we 
discuss and conclude our results. 



II. GALACTIC CLUSTERS IN THE 
SCALAR-TENSOR VERSION OF f(R) GRAVITY 

Astronomical observations have shown that galaxies 
tend to concentrate in larger structures, called clusters of 
galaxies. The total mass of galaxy clusters ranges from 
10 13 Mq for groups up to a few 10 15 M & for very rich sys- 
tems, and the cluster morphology is usually dominated 
by a regular centrally peaked main component [H, |36| . 
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Since clusters are "dark matter" dominated objects their 
formation and evolution is driven by gravity. In this con- 
text, the mass function of the clusters is determined by 
the initial conditions of the mass distribution set in the 
early universe, and the evolution of the large scale mat- 
ter distribution on scales comparable to the size of the 
clusters is linear. The overall process of the gravitational 
growth of the density fluctuations and the development 
of gravitational instabilities leading to cluster formation 
has been extensively studied by using both analytical and 
numerical methods [37| . 

In order to derive the generalization of the relativis- 
tic virial theorem for galaxy clusters in the f(R) gravity 
models we need, as a first step, to obtain the gravita- 
tional field equations for a static and spherically sym- 
metric system, with the matter energy-momentum tensor 
corresponding to a system of identical, collisionless point 
particles in random motion. To derive the basic field 
equations we will use the scalar-tensor representation of 
f(R) gravity, which allows a clear physical interpretation 
of the model. 



A. Scalar-tensor representation for f(R) gravity 

The action for f(R) modified theories of gravity takes 
the form 



where 



S = 



1 



16ttG 



f(R) + L Tl 



(1) 



where f(R) is an arbitrary analytical function of the Ricci 
scalar R and L m is the Lagrangian density corresponding 
to matter. 

Varying the action with respect to the metric g^ yields 
the following field equations 

FR»u ~ \f9^ - (V M V„ - 9tiV U) F = BnGT^, (2) 

where we have denoted F — df /dR and T M „ is the stan- 
dard minimally coupled matter energy-momentum ten- 
sor. Note that the covariant derivative of the field equa- 
tions and of the matter energy-momentum tensor van- 
ishes for all f(R) by means of generalized Bianchi iden- 
tities [H, H|| . Contracting the field equations gives the 
useful relation 



3DF + FR — 2/ = 8nGT, 



(3) 



from which one verifies that the Ricci scalar is now a fully 
dynamical degree of freedom, and where T is the trace of 
the energy-momentum tensor. 

By introducing the Legendre transformation {R, /} — > 
{4>,V} defined as 

<t> = F(R), V(4>) = R^)F- f(R(<f>)), (4) 

the field equations of f(R) gravity can be reformulated 
as 39] 



1 G 
Riiv — -^g^vR = Sk—T^ 



(5) 



~V (</>) g„u + \ (V M V„ - g^D) 4>. 



Using these variables, Eq. takes the form 
3D(j> + 2V(<j)) - <j>^- = 8nGT. 



(6) 



(7) 



In this representation the field equations of f(R) grav- 
ity can be derived from a Brans-Dicke type action with 
parameter u = 0, given by 



S = 



1 



16ttG 



W - v{4>) 



-g d 4 x. (8) 



The only requirement for the f(R) model equations 
to be expressed in the form of a Brans-Dicke theory is 
that F(R) be invertible, that is, R(F) exists ^20]. This 
condition is necessary for the construction of V(<f>). 

The modification of the standard Einstein-Hubert ac- 
tion leads to the appearance in the field equations of an 
effective gravitational constant G e ff = G/<f>, which is a 
function of the curvature. Secondly a new source term 
for the gravitational field, given by the tensor 0^ v , is also 
induced. The tensor 9^ v is determined by the trace of the 
energy-momentum tensor via Eq. ([7]) , which thus acts as 
an independent physical parameter determining the met- 
ric of the space-time. 



B. Field equations for a system of identical and 
collisionless point particles 

Let us assume a time-oriented Lorentzian four- 
dimensional space-time manifold M.. Consider an iso- 
lated and spherically symmetric cluster described by a 
static and spherically symmetric metric 



ds 2 



i v ^dt 2 + e^dr 2 + r 2 (d0 2 + sm 2 0dcp 2 ) . (9) 



The galaxies, which are treated as identical and col- 
lisionless point particles, are described by a distribution 
function fg. The latter distribution function obeys the 
general relativistic Boltzmann equation, which shall be 
treated in detail below. 

In terms of the distribution function the energy- 
momentum tensor of the matter can be written as 14011 



T, 



f B m u^Uy du, 



(10) 



where m is the mass of the particle (galaxy), is the 
four- velocity of the galaxy and du = du r duedu v / Ut is the 
invariant volume element of the velocity space. There- 
fore, the energy-momentum tensor of the matter in a 
cluster of galaxies can be represented in terms of an effec- 
tive density p e s and of an effective anisotropic pressure, 
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with radial p e r 2 and tangential p^g components, given 
by 



PeS 



^-p(4) = p{ul) 



(11) 



where, at each point, (u 2 ) is the average value of u 2 , etc, 
and p is the mass density [32j . 

By using this form of the energy- momentum tensor, 
the gravitational field equations describing a cluster of 
galaxies, in f(R) gravity, take the form 



III. THE VIRIAL THEOREM IN f(R) 
GENERALIZED GRAVITY MODELS 

In order to derive the virial theorem for galaxy clus- 
ters, which are described by the distribution function h, 
we have to first write down the Boltzmann equation gov- 
erning the evolution of the distribution function. This 
equation can then be integrated over the velocity space, 
to yield an equation which, used in conjunction with the 
gravitational field equations, provides the required gen- 
eralization of the virial theorem. 



- x I 1 j 1 g ^ / 2 \ 

r 2 r J r 2 <h^ 



1 

2d> 



V(<f>) 



1 



(V t V -□)</>, (12) 



-e~ A I i' 



a) 4, (13) 



,1- 



X' 



G 



1 



tor- p (*») + — V{4>)--{V e V e -U)cp (14) 



(V v V v -D)0. (15) 



A useful relationship is obtained by adding the grav- 
itational field equations Eqs. (fl"2]) - (|15p . from which we 
obtain the following equation 



v 
r 



v'\' 



+ V(0) + i(2V t V* + □)</>, (16) 

where (u 2 ) = (u 2 ) + (u 2 r ) + (u 2 ) + {u%). 

Since we are interested in astrophysical applications at 
the extra-galactic level, we may assume that the devia- 
tions from standard general relativity (corresponding to 
the background value = 1) are small. Therefore we may 
represent 4> as 4> — 1 + e <?'(-R)j where e is a small quantity, 
and g'{R) describes the modifications of the geometry 
due to the presence of the tensor 0„„ ,20]. Consequently 
I /(f> ~ 1 — eg'(R), and Eq. (|16[) can be rewritten as 



— A 



V 

r 



v'X 1 



AnGp{u 2 ) + AirGpcf,, 
(17) 



where 
47rGp0 



-4irGep(u 2 )g'{R) 



1 



-V ((/))+ - (2V t V 



□ 



(18) 



A. The relativistic Boltzmann equation 

Consider a tangent bundle T(A4), which is a real vec- 
tor bundle whose fibers at a point x £ M. is given by 
the tangent space T x (M). In the space-time M the in- 
stantaneous state of a particle with mass mo is given by 
a four- momentum p S T x (M) at an event x €E M.. The 
one-particle phas e sp ace Pphase is a subset of the tangent 
bundle given by [401 ] 

pphase := {(x,p) | x £ M,p e T x (M) ,p 2 = -to 2 ,} . 

(19) 

A state of a multi-particle system is described by a 
continuous, non-negative distribution function /s(x,p), 
defined on Pphase, and which gives the number dN of 
the particles of the system which cross a certain space- 
like volume dV at x, and whose 4-momentap lie within a 
corresponding three-surface element dp in the momentum 
space. The mean value of fs gives the average number of 
occupied particle states {x,p). Macroscopic, observable 
quantities can be defined as moments of fs [401 ] . 

Let {x Q }, a = 0, 1, 2, 3 be local coordinates in an open 
set U C A4. The coordinates are chosen so that dt is 
time-like future directed and d a , a — 1,2,3, are space- 
like. Then {d/dx a } is the corresponding natural basis 
for tangent vectors. We express each tangent vector p in 
U in terms of this basis as p — p a d/dx a and define a sys- 
tem of local coordinates A — 0, . . . , 7 in Tjj(M) as 
z a = x a , z a+A = p a . This defines a natural basis in the 
tangent space given by {d/dz A } = {d/dx a , d/dp a } (40T |. 

A vertical vector field over T (M) is given by tt = 
p a d/dp a . The geodesic flow field a, which can be con- 
structed over the tangent bundle, is defined as a = 

where are the 



p a d/dx a - p a P T^d/dp 3 = p a D a 
connection coefficients. Physically, a describes the phase 
flow for a stream of particles whose motion through 
space-time is geodesic [40j. 

Therefore, the transport equation for the propagation 
of a particle in a curved arbitrary Riemannian space-time 
is given by the Boltzmann equation [40( 



b=l+eg'(R) 



dx a 



p a p Y 



d_ 



fs = 0. 



(20) 
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For many applications it is convenient to introduce 
an appropriate orthonormal frame or tetrad e.®(x), a — 
0,1,2,3, which varies smoothly over some coordinates 
neighborhood U and satisfies the condition g^e^e^, = 
rj ab for all x € U [H, H(| • Any tangent vector p^ at x 
can be expressed as p^ = p a e£, which defines the tetrad 
components p a . 

In the case of the spherically symmetric line element 
given by Eq. (J9j) we introduce the following frame of or- 
thonormal vectors [33, Sol : 



o = v /2 8 o 



rS 2 

' "(I! 



e 1 = e A/2 <5 1 
e 3 , = r sin 95 3 , 



(21) 



Let u M be the four-velocity of a typical galaxy, sat- 
isfying the condition u^u^ = — 1, with tetrad compo- 
nents u a — u^e"^. The relativistic Boltzmann equation in 
tetrad components is 

a dx^ bc du a 



0. 



(22) 



where the distribution function fg = Jr(x^, u a ) and 
Ibc = e p;f e fc e c are the Ricci rotation coefficients (32l.l40l|. 
By assuming that the only coordinate dependence of the 
distribution function is upon the radial coordinate r, 
Eq. (JUD becomes ^ 



dr 



1 2 dv 
1 

wi 



r 

dfs 



U-2 



du 2 



dfs 
du\ 

dfg 
du 3 



u 3 



-e x/2 u 3 cote f 



df B U6 9f B \ 
du 3 du 2 J 



= 0. (23) 



B. The generalized virial theorem in f(R) modified 
theories of gravity 



The spherically symmetric nature of the problem re- 
quires that the coefficient of cot#, in Eq. (|23|) . be zero, 
which implies that the distribution function fg is only a 
function of r, u% and u\ + u\. Now, multiplying Eq. (j2"3")l 
by mu r du, then integrating over the velocity space, and 
by assuming that fg vanishes sufficiently rapidly as the 
velocities tend to ±oo, we obtain 

-p[{ul) + {ul)-2{u\)] =0. (24) 

The following step consists in multiplying Eq. (j24|) by 
47rr 2 , and integrating over the cluster provides the fol- 
lowing relationship [32j 



4irp [(uf) + (ul) + (ul)} r 2 dr 



dv 



\itr 3 p [(u 2 ) + (u 2 )] —dr = 0. (25) 



It is convenient to introduce some approximations that 
apply to test particles in stable circular motion around 
galaxies, and to the galactic clusters. First of all, we 
assume that v and A are slowly varying (i.e. v' and A' 
small), so that in Eq. (fl~7|) the quadratic terms can be 
neglected. Secondly, we assume that the galaxies have 
velocities much smaller than the velocity of the light, so 
that (u\) « (ul) « (ul) < (ul) w 1. Thus, Eqs. (JT7J) 
and (|25p become 

i£( r2 £) =47rGp + 47rG ^' (26) 



and 



2K--[ A^p^-dr = 0, 
2 J dr 



(27) 



respectively, where 
r-rt 



K 



2np [{u 2 ) + (u 2 ) + (ul)] r 2 dr, (28) 



is the total kinetic energy of the galaxies. The total mass 
of the system is given by M = dM(r) = J^ 1 Airpr 2 dr. 
Note that the main contribution to M is due to the bary- 
onic mass of the intra-cluster gas and of the stars, but 
other particles, such as massive neutrinos, may also con- 
tribute significantly to M. 

Now, multiplying Eq. (|26|) by r 2 and integrating from 
to r we obtain 



GM{r)= l -r 2 ^--GM^(r) 



where we have defined 



M (r) =4^ / P4 ,(r')r u dr' 
Jo 



(29) 



(30) 



We may call this useful quantity as the geometric mass 
of the cluster. By multiplying Eq. (|2"9"j) with dM(r), fol- 
lowing an integration from to R, and by defining 



GM(r) 



dM(r), 



and 



Jo 



r 



we obtain the following relationship 



n = n 



4> 



1 

2 Jo 



'a 3 dv j 

iirr p— dr, 
or 



(31) 



(32) 



(33) 



where the quantity Q is the usual gravitational potential 
energy of the system. 

Finally, with the use of Eq. (|27|) , we obtain the gener- 
alization of the virial theorem, in f(R) modified theories 
of gravity, which takes the form 



2K + n - fls = 0. 



(34) 
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Note that the generalized virial theorem, given by 
Eq. (|34p , can be written in an alternative form if we in- 
troduce the radii Ry and defined by 



and 



R\ 



R<i 



M 



R M{r) 



dM(r), 



M. 



r 



(35) 



(36) 



respectively. We denote i?^ as the geometric radius of 
the cluster of galaxies. Thus, the quantities and 
are finally given by 



[2 



and 



GM 
Rv 



GM, 



r2 



R. 



■<t> 



(37) 



(38) 



respectively. 

The virial mass My is defined as f4li 



2K 



GMM\ 
Rv 



(39) 



After substitution into the virial theorem, given by 
Eq. Ij3"3]). we obtain 



My _ M\R V 
M M 2 R th ' 



(40) 



If My/M > 3, a condition which is valid for most of 
the observed galactic clusters, then Eq. |40|) provides the 
virial mass in f(R) gravity, which can be approximated 
by 



M v 



MlRv 
M R 6 ' 



(41) 



The virial theorem can also be derived in the frame- 
work of Newtonian gravity, and it can be formulated as 



d 2 i. 



jk 



dt 2 



= AK jk + 2W ik 



-A/ 



jk 



(42) 



where Ijk is the moment of inertia tensor of the gravitat- 
ing system, Kjk is the kinetic energy tensor, Wjk is the 
gravitational potential energy tensor and A is the cosmo- 
logical constant. If an external force with potential <& cxt 
acts on the system, a new term of the form 

V jk =--Jp [x h -^- + j dV, (43) 

must be added to the right-hand side of Eq. (|4"2l . In 
gravitational equilibrium, when d 2 Ijk/dt 2 = 0, the virial 



theorem gives 2K — \W\ — 0, where K and W are the 
traces of Kjk and Wjk, respectively, and we have as- 
sumed A = 0. For \W\ = (3/5)GM v /R v , correspond- 
ing to a constant density matter distribution, we obtain 
that K = (3/10)GM 2 /Ry. However, in the "geometric 
mass" approach to the dark matter problem, in which 
one has to make a clear distinction between matter, in 
the usual sense, and geometric effects, only the baryonic 
matter has kinetic energy, which requires a modification 
of the definition of the virial mass. Moreover, the geo- 
metric effects, specific to f{R) gravity models, cannot be 
systematically considered in a purely Newtonian deriva- 
tion. 

In the observational context, the virial mass My is 
determined from the study of the velocity dispersion a 2 
of the stars and of the galaxies in the clusters. Now, 
according to our interpretation, most of the mass in a 
cluster with mass M tot should be in the form of the ge- 
ometric mass M<f„ so that M^ « M tot . A possibility 
of detecting the presence of the geometric mass and of 
the astrophysical effects of the /(i?) extensions of gen- 
eral relativity is through gravitational lensing, which can 
provide direct evidence of the mass distribution and of 
the gravitational effects even at distances extending far 
beyond of the virial radius of the galaxy cluster. 



IV. ASTROPHYSICAL APPLICATIONS 



It is important to emphasize that astrophysical ob- 
servations together with cosmological simulations have 
shown that the virialized part of the cluster corresponds 
roughly to a fixed density contrast S ~ 200 as compared 
to the critical density of the universe, p c (z), at the con- 
sidered redshift, so that pv — 3 My /4irR v = Spdz), 
where pv is the virial density, My and Ry are the virial 
mass and radius; p c is given by p c (z) = h 2 {z)3HQ /8ttG, 
where h(z) is the Hubble parameter normalized to its lo- 
cal value, i.e., h 2 iz) — fl m (1 + z) + where fl m is 
the mass density parameter and J7a is the dark energy 
density parameter (35j . 

Now, once the integrated mass as a function of the ra- 
dius is determined for galaxy clusters, a physically mean- 
ingful fiducial radius for the mass measurement has to be 
defined. The radii commonly used are either r 2 oo or r 50 o- 
These radii lie within the radii of the mean gravitational 
mass density of the matter (ptot) = 200p c or 500p c . A 
pragmatic approach to the virial mass is to use r2oo as 
the outer boundary [361 ]. The numerical values of the 
radius r2oo are in the range r2oo = 0.85 Mpc (for the 
cluster NGC 4636) and r 200 = 4.49 Mpc (for the cluster 
A2163), so that a typical value for r2oo is approximately 
2 Mpc. The masses corresponding to r2oo and r^oo are 
denoted by M200 and Af5oo, respectively, and it is usually 
assumed that My = M200 and Ry = r 2 oo [III- 
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A. Geometric mass and geometric radius from 
galactic cluster observations 

Most of the baryonic mass in the clusters of galaxies is 
in the form of the intra-cluster gas. The gas mass density 
p g distribution can be fitted with the observational data 
by using the following expression for the radial baryonic 
mass (gas) distribution [36| 



Pg( r ) =P0 1 + 



-3/3/2 



(44) 



where r c is the core radius, and po and (3 are (cluster- 
dependent) constants. 

It is usually assumed that the observed X-ray emis- 
sion from the hot, ionized intra-cluster gas is in isother- 
mal equilibrium. Therefore, one may assume that the 
pressure P g of the gas satisfies the equation of state 
Pg = {k B T g / pm p )p gi where k B is Boltzmann's constant, 
T g is the gas temperature, p 0.61 is the mean atomic 
weight of the particles in the cluster gas, and m p is the 
proton mass [3a]. Thus, with the use of the Jeans equa- 
tion [lj we obtain the total mass distribution as [34], [3(| 



M tot {r) 



knT, 



9 r 2 



dr 



In / 



(45) 



pm p G 

Now, taking into account the density profile of the gas 
given by Eq. (14~4"|) , we obtain for the total mass profile 
inside the cluster the following relation [361 ] 



M tot (r) 



3k B f3T g 



(46) 



pm p G r 2 , , 

According to the modified f(R) gravity model, the to- 
tal mass of the cluster consists of the sum of the baryonic 
mass (mainly the intra-cluster gas), and the geometric 
mass, so that M tot (r) = 4ir JJ* [p g + p^r^dr. Hence it 
follows that M to t (f) satisfies the following mass continu- 
ity equation, 

dM tot (r) 



dr 



Aitr pg (r) + A-kt p^ (r) . 



(47) 



Since the gas density and the total mass profile inside 
the cluster are given by Eqs. ([4*4")) and (|4*6"|) , respectively, 
we can immediately obtain the expression of the geomet- 
ric density term inside the cluster as 

y.2 



4ir P4> (r) 



3k B PT g (r 2 + 3rg) 
pm p (r 2 + r 2 ) 2 
In the limit r > r c wc obtain for 
3k B f3T g 



AnGpo 



. . /2 - (48) 
(l + r 2 /r c 2 ) 3/3/2 

30 the simple relation 



47rp (r) 



prrip 



4rtGp r, 



3^2-3/3 



(49) 



The geometric mass can be obtained generally as 



GM,p(r) — 47r / r 2 p^{r)dr 
Jo 

- AirGpc 



Sk B f3T q 



pm p 1 + r 2 /r 2 
r 2 dr 

(l + r 2/ r 2)3/3/2' 



(50) 



and in the limit r » r c , may be approximated as 



GM$ (r) 



3k B /3T g 4irG Po r 



pm p 



3/3„2-3/3 



3(1-/3) 



(51) 



One may assume that the contribution of the gas den- 
sity and mass to the geometric density and geometric 
mass, respectively, can be neglected. The latter approx- 
imations are very well supported by astrophysical obser- 
vations, which show that the gas represents only a small 
fraction of the total mass [351 l36t [37l |42| . Therefore, we 
obtain 



AirGp^r) 



3k B [3T, 



pm p 



and 



(52) 



(53) 



respectively. 

By using the obtained mass profile one can obtain the 
form of the metric tensor component e 1 ' inside the clus- 
ter. From Eq. (|2"6"1) it follows that v' must satisfy the 
condition r 2 v' = 2M to tir) f» 2M§ (r), which gives 



(54) 



where C v is an arbitrary integration constant, and we 
defined 



3k B pT g 
pGm p ' 



(55) 



for notational simplicity. As for the metric coefficient 
e~ A , we can assume that in the first approximation it is 
given by its standard general relativistic form, e~ A k1- 
2GM tot /r, and is given by e _A »1- %k B (3T g / pm p = 1 — 
2Gs. However, this approximation may not be necessary 
correct inside the clusters of galaxies. 

Hence, in modified f(R) gravity models the metric in- 
side a galactic cluster can be directly obtained from as- 
trophysical observations. 

One may also estimate an upper bound for the cut- 
off of the geometric mass. The idea is to consider the 
point at which the decaying density profile of the geo- 
metric density associated to the galaxy cluster becomes 
smaller than the average energy density of the Uni- 
verse. Let the value of the coordinate radius at the 
point where the two densities are equal to be R^ r ^ ■ 

Then at this point p^{R^ r) ) = p univ , where p univ is 
the mean energy density of the universe. By assuming 
Pumv = Pc = 3ff 2 /87rG = 4.6975 x 10- 30 /i| g/c m - 3 , 
where H = 50/iso km/Mpc/s (36|, we obtain 



R 



(cr) 



/ 3k B {3T g 



1/2 



\pm p Gp, 
= 9L33 ^5keVj 



k B Tg\ 



h^Mpc. (56) 
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The total geometric mass corresponding to this value 



is 



4.83 x 10 16 /? 3 / 2 (jg^ j ^M . (57) 



fk B T g \ 3/2 



This value of the mass is consistent with the observations 
of the mass distribution in the clusters of galaxies. How- 
ever, according to f(R) modified theories of gravity, we 
predict that the geometric mass and its effects extends 
beyond the virial radius of the clusters, which is of the 
order of only a few Mpc. 



B. Radial velocity dispersion in galactic clusters 

The virial mass can also be expressed in terms of the 
characteristic velocity dispersion 01 as [42| 



My = -a\R v . 



(58) 



By assuming that the velocity distribution in the cluster 
is isotropic, we have (u 2 ) = {u\) + (u 2 ) + (uf) = 3(u 2 ) = 
3er 2 , where er 2 is the radial velocity dispersion, o\ and 
oy are related by 3er 2 = <r 2 . 

In order to derive the radial velocity dispersion rela- 
tion for clusters of galaxies in f(R) gravity we start from 
Eq. Taking into account that the velocity distribu- 

tion is isotropic, we obtain 



1 dv 



(59) 



One may assume that inside the cluster the condi- 
1, to a first order of approximation, and in 



tion e 



-A 



the limit of small velocities the modified field equation, 
Eq. (fTT|) , may be integrated to yield 



As an example of the application of Eq. (|62|). consider 
the case in which the density p of the normal matter 
inside the cluster has a power law distribution, so that 



p(r) = p r 7 , 



(63) 



1,3 positive constants. The cor- 
matter mass profile is M(r) — 
The geometric mass GM^ sa 

is 



with po and 7 ^ 
responding normal 
4 W 3 -V(3-7)- 

(SksftTg/ prrip) r = qor, where qo = (3kB(3T g / pm 
linearly proportional to r, as has been shown in the pre- 
vious Section. 

Therefore, for 7^1,3, we obtain the following solution 



qo 



a r (r) = — 



2nGp 



r 2- 7 



C 1 d 



7 (7-I) (3 -7) 



7 + 1 r 



Po 



For the specific case of 7 = 1, we find 
and for 7 = 3, we have 



C C\ 

q + 2TrGp r In r H r, 

2r po 



2 1 \ 1o ri 



In / 



1 



CI 

4 r 



— '/ 



(64) 



(65) 



(66) 



The observed data for the velocity dispersion in clus- 
ters of galaxies are usually analyzed by assuming the sim- 
ple form <7 2 (r) = B/(r + b) for the radial velocity disper- 
sion, with B and b constants. As for the density of the 
galaxies in the clusters the relation p(r) — A/r (r + a) 2 , 
with A and a constants, is used. The data arc then fitted 
with these functions by using a non-linear fitting proce- 
dure [42j |. For r a, p(r) A/r, while for r a, 
p(r) behaves like p{r) w A/r 3 . Therefore the comparison 
of the observed velocity dispersion profiles of the galaxy 
clusters and the velocity dispersion profiles predicted by 
f(R) modified theories of gravity may give a powerful 
method to discriminate between the different theoretical 
models. 



rV 



2GM 4> (r) + 2GM(r) + 2C, 



(60) 



where C is an arbitrary constant of integration. Since 
from Eq. ([59^1 we have v' = -(2/ p)d (pa 2 ) /dr, it follows 
that the radial velocity dispersion of the galactic clusters, 
in f(R) gravity, satisfies the differential equation 



dr 



\po r ) = 



GM^r) GM(r) C 
~2 P( r ) —P( r ) - r*P( r )> 



(61) 



and which in turn provides the following general solution 



a 2 (r) = - 



GM^r') GM(r') . , 



-Pin 



p(r 

rfr' + — , (62) 
P 



where C\ is an integration constant. 



THE LAGRANGIAN OF THE f(R) GRAVITY 
MODEL 



The virial theorem in /(i?) gravity models, which leads 
to the possibility of obtaining the metric tensor compo- 
nents inside the cluster, also opens the possibility of di- 
rectly obtaining the Lagrangian f(R) of the theory from 
astrophysical observations. To obtain the Lagrangian we 
start from the field equations (fL^ - (fT5")) . which in the stan- 
dard representation can be written as 



F" 



1 



(1/ + A') F' 



A' 



-F = 0, 



(67) 



," + ,' 2 + A') -l(l-e^) 



-2 F - 
F 



(68) 
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f = Fe 



and 



-A 



v"-\{v' + \')v'--\' + + 

2 r \ r J 



~F 



(69) 



F I F 



2 r 



T> 171)1 



respectively, where we have neglected the contribution of 
the baryonic matter (intra-clustcr gas, stars, etc). By 
using for exp(t^) the expression given by Eq. (|54[) , and 
by assuming that inside the cluster exp(A) =constant, 
Eq. (l6~7| gives a second order linear differential equation 
for F(r), 



r 2 F" - srF' - 2sF = 0, 
with the general solution 

F(r) = C ir q - +C 2 r q+ , 



(71) 



(72) 



where C\ and Ci are arbitrary integration constants, 
and q± = (l + s ± Vl + 10s + s 2 ) /2. By consid- 
ering for F only the monotonically decreasing solu- 
tion we obtain F(r) — C\r~ q , where q = — g_ = 
[Vl + 10s + s 2 - (1 + s)] j2 > 0. With this form of the 
function F(r), Eq. (|6"8")) fixes exp(A) as 



l + s(2-s)-q(q + 2) 



= 1 - 8s - 3s 2 + 2s\/l + 10s + s 2 



(73) 



Eq. (jnij) gives the radial coordinate dependence of the 
Lagrangian / as 



f = 2de- x 4(1 + s) -2(2 + s)VT + T0s+7 2 

while from Eq. (|70[) we obtain for R the expression 

„ (1 + 5s) Vl + 10s + s 2 - 1 - 34s - 9s 2 _ A 
R = o e 



9-2 



(74) 



(75) 



Finally, we obtain / as a function of R in the form 

f(R) = f R 1+q/2 , (76) 

where 



fo = 2Ci x 



4(1 + s) - 2(2 + s)Vl + 10s + s 2 



[(1 + 5s) VI + 10s + s 2 - 1 - 34s - 9s 2 ] 



1+9/2 



(77) 



Therefore, once the main physical parameters of the 
gas in the cluster, such as the gas temperature T g or the 
gas density profile, described by the parameter (3, are 
known, the action of the modified gravity model can be 
completely obtained from observations. 

The weak field limit of the f(R) generalized gravity 
models has been discussed recently, for star-like objects, 



in [18| and [l£|, respectively. By assuming that f(R) 
is an analytical function at the constant curvature Rq, 
that m^r <§; 1, where is the effective mass of the 
scalar degree of freedom of the theory, and that the 
fluid is pressureless, the post-Newtonian potentials \& (r) 
and $ (r) are obtained for a metric of the form ds 2 - 



[1 - 2* (r)] dt 2 + [1 + 2$(r)] dr 2 



(d9 2 



] 9d V 2 ). 



One may then find the behavior of "J (r) and $(r) outside 
the star. This analysis leads to a value of 7 = 1/2 for 
the Post-Newtonian parameter 7, which from Solar Sys- 
tem observations is known to have a value of 7 = 1. An 
analysis of a Lagrangian of the form given by Eq. (|76[) . 
denoted in [l8[ as f(R) = (R/a) 1+s has also been con- 
sidered, and the conclusion is that ". . .this analysis is 
incapable of determining whether f(R) = R 1+s gravity 
with 6 7^ 1 conflicts with Solar System tests" [l8[ . There- 
fore, once the main physical parameters of the gas in the 
cluster, like the gas temperature T g , or the gas density 
profile, described by the parameter /?, are known, q can 
be calculated, the action of the modified gravity model 
can be completely obtained from observations, and the 
viability/non- viability of the model can be directly tested 
by using cluster of galaxy data, which may offer an al- 
ternative to the Solar System tests. 

Another problem facing the f(R) gravity models is the 
problem of the stability. On a time scale of r w 10~ 26 
s a "fatal instability" develops when f"(R) < [l!|. 
For the Lagrangian given by Eq. (fTB")) we have f"{R) = 
fo(q/2)(q/2 + ^Rt' 2 ' 1 > 0. Therefore this type of in- 
stability does not develop in the present model. 



VI. SUMMARY AND DISCUSSIONS 

Cosmology has entered a 'golden age' in which the 
rapid development of increasingly high-precision data has 
turned it from a more theoretically driven to an obser- 
vationally based science. Nevertheless, modern astro- 
physics and cosmology are plagued with two severe dif- 
ficulties, known as the dark energy and the dark mat- 
ter problems. One promising approach of theoretical 
research to improve our understanding of these issues 
is modified gravity. In particular, f(R) modified theo- 
ries of gravity challenge the need for dark matter and 
for dark energy. Some of the models seem to account 
for the late time acceleration of the universe, and vi- 
able models seem to exist. They simultaneously account 
for the four distinct cosmological phases, inflation, the 
radiation-dominated and matter-dominated epochs, and 
the late-time accelerated expansion, respectively. More- 
over, these models also seem to be consistent with cos- 
mological structure formation observations. 

In an astrophysical context, galactic dynamics of mas- 
sive test particles may also be understood without the 
need for dark matter in the framework of f{R) modified 
theories of gravity. It is therefore of paramount inter- 
est to derive a generalized version of the virial theorem, 
which plays an important role in many areas of astro- 
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physics, in the context of the modified theories of gravity. 
The virial theorem is an extremely useful tool in stabil- 
ity issues, and in the deduction of the mean density of 
astrophysical objects such as galaxies, clusters and super 
clusters, and consequently in the prediction of the total 
mass of the clusters of galaxies. 

In this work, we have analyzed the 'dark matter' prob- 
lem by considering a generalized version of the virial the- 
orem in the framework of f(R) modified theories of grav- 
ity. The virial theorem was obtained by using a method 
based on the collisionless Boltzmann equation. The ad- 
ditional geometric terms present in the modified gravi- 
tational field equations provide an effective contribution 
to the gravitational energy, which at the galactic/extra- 
galactic level acts as an effective mass, playing the role 
of the 'dark matter'. The total virial mass of the galac- 
tic clusters is mainly determined by the effective mass 
associated to the new geometrical term, the geometrical 
mass. It is important that the latter term may account 
for the well-known virial theorem mass discrepancy in 
clusters of galaxies. 

In the present model there is also a strict proportional- 
ity between the virial mass of the cluster and its baryonic 
mass, a relation which can also be tested obscrvationally. 
Since galaxy clusters are "dark" matter dominated ob- 
jects, the main contribution to their mass comes from 
the geometric mass M^, so that with a very good ap- 
proximation we have as My rs M to t- Therefore the 
virial theorem gives immediately the following mass scal- 
ing relation 

M v « M^. (78) 
Rv 

This equation shows that the virial mass is proportional 
to the baryonic (normal) mass of the cluster, and that 
the ratio of the total mass and of baryonic mass is de- 
termined by a purely geometric quantity, the geometric 
radius R^. Hence the geometric radius of the cluster 
can be determined from observations, once the virial and 
baryonic masses and the virial radius, respectively, are 
known. 

By assuming that R$ ~ Rl , , we obtain the following 
relation between the virial and the baryonic mass of the 
cluster 

Mv , 91 .33V?(|^) 1/2 W^. (79) 



For a cluster with gas temperature T g — 5 x 10 7 K, j3 = 
1/2 and Rv — 2 Mpc we obtain My as 32M, a relation 
which is consistent with the astronomical observations 




One of the important predictions of the f(R) "dark 
matter" model is that the geometric mass and its effects 
extends beyond the virial radii of the clusters. Gener- 
ally, the virial mass Mv is obtained from the observa- 
tional study of the velocity dispersions of the stars in the 
cluster. Due to the observational uncertainties it can- 
not give a reliable estimation of the numerical value of 
the total mass M g + M^ in the cluster. However, a much 
more powerful method for the study of the total mass dis- 
tribution in clusters is the gravitational lensing of light, 
which may provide direct evidence for the gravitational 
effects at large distances from the cluster. The presence 
of f{R) modified gravity effects at large distances from 
the cluster leads to significantly different lensing observa- 
tional signatures, as compared to the standard relativis- 
tic/dark matter model case. The bending angle in the 
f(R) models could be larger than that predicted by the 
dark matter models. Therefore, the observational study 
of the gravitational lensing could discriminate between 
the different dynamical laws proposed to model the mo- 
tion of particles at the clusters of galaxy level, and the 
standard dark matter models. 

In conclusion, the generalized virial theorem in f(R) 
gravity might be an efficient tool in observationally test- 
ing the viability of this class of generalized gravity mod- 
els. 
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